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The connection between the thermodynamics of charged finite nuclear systems and the 
asymptotically measured partitions is presented. Some open questions, concerning in par- 
ticular equilibrium partitions are discussed. We show a detailed comparison of the decay 
patterns in Au + C, Cu, Au central collisions and in Au quasi-projectile events. Obser- 
vation of abnormally large fluctuations in carefully selected samples of data is reported 
as an indication of a first order phase transition (negative heat capacity) in the nuclear 
equation of state. 

><j : 1. INTRODUCTION 



In the last 20 years and especially in the most recent ones, considerable progress has 
been achieved both theoretically and experimentally in the investigation of nuclear reac- 
tion dynamics and thermodynamics in the Fermi energy regime [1]. 

From a theoretical point of view, strong efforts have been devoted to the understand- 
ing of the nuclear equation of state (EoS) either with transport theories or statistical 
approaches. 

From an experimental point of view new generation 47r detectors have been developed 
and are now operating at different accelerator facilities (Dubna, GANIL, GSI, LNL, LNS, 
MSU, Texas A-M). They are producing a huge amount of exclusive data and new kind 
of analyses. Very rich information has already been extracted from experimental stud- 
ies on intermediate energy heavy ion collisions. Experiments have shown that the final 
state can be constrained to select the dynamics of the collision and isolate events that 
populate states closely compatible with equilibrium. Several investigations have demon- 
strated that excited nuclear systems produced in such collisions undergo in a short time 
scale (100 fm/c) bulk multifragmentation characterized by final states containing several 
Intermediate Mass Fragments (1MF,Z > 3). 

A considerable progress has been accomplished on the theoretical as well as on the ex- 
perimental side in order to define and collect a converging ensemble of signals connecting 
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multifragmentation to the nuclear liquid-gas phase transition [1] and locating it in the 
phase diagram of nuclear matter [2]. The opening of the high fragment multiplicity chan- 
nel, the onset of collective expansion, the flattening of the caloric curves, the fossile signal 
of spinodal decomposition, a negative branch of the heat capacity, the bimodal distribu- 
tion of exclusive observables and the finite size and Fisher law scalings have been observed 
and tentatively related to the Equation of State of the nuclear matter. The possibility 
of new radioactive beams (RIB) facilities is now prompting exciting theoretical advances 
of the isospin aspects of the EoS, like the density dependence of the symmetry energy 
and the modifications to the spinodal instability [1]. Experimentally, first studies on the 
dependence of the isoscaling parameters on the isospin of the decaying system [3] already 
started exploiting stable beams. These works seem to indicate an isospin distillation in 
asymmetric systems. 

All these signals can be considered as circumstantial evidences of a phase transition, but 
some of them are still controversial and need to be further experimentally investigated 
before the phase transition can be definitely assessed. This work is a contribution to 
this aim, but only next generation experiments will allow to reach a comprehensive and 
detailed understanding of the phase transition and the nuclear EoS. 

2. THERMODYNAMICS OF NUCLEAR SYSTEMS 

2.1. EQUILIBRIUM PARTITIONS 

In order to perform thermodynamical analyses, one has to collect a data sample which 
corresponds as closely as possible to an homogeneous population of the phase space. 
Data must be selected such as to isolate a portion of the cross section where the entire 
system (or the quasi-projectile) properties keep a negligeable memory of the entrance 
channel dynamics. This can be experimentally verified checking that for a given source 
the fragmentation pattern is determined by the size, charge, energy and average freeze- 
out volume solely, independent of the way the source has been formed, e.g. different 
impact parameters. In this case the thermodynamics we can access is a microcanonical 
thermodynamics with energy, number of protons and neutrons, and average volume as 
state variables. Indeed the excitation energy can be measured on an event-by-event basis 
by calorimetric techniques. For any shape of the excitation energy distribution the events 
can thus be sorted in constant energy bins, i.e. in microcanonical ensembles. 

In the following we present thermodynamical studies performed on quasi-projectile 
events from peripheral 35 A MeV Au + Au collisions [4] and central events from 25 A MeV 
Au + C, 25 and 35 A MeV Au + Cu and 35 A MeV Au + Au collisions [5], measured at the 
K1200-NSCL Cyclotron of the Michigan State University with the MULTICS-MINIBALL 
apparatus. 

Single source almost complete events have been selected with a constant value for the 
collected charge 1 in each energy bin [4,5]. The possible pollution from other sources has 
been minimized for central collisions through a shape analysis and in the case of Au 
quasi-projectile by substituting the backward light particle emission by the symmetric of 
the forward emission in the quasi-projectile reference frame. The observed event isotropy 
(Fig. 1) indicates that the directed flow component coming from a memory of the entrance 

1 90% of the total charge for central collisions and 90% of the projectile charge for peripheral events 
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Figure 1. Left: Angular distribution of the fragments (Z > 3), but the heaviest one, in the Au 
quasi-projectile reference frame in different intervals of the excitation energy. Right: Angular 
distribution in the center of mass reference frame for central Au + Cu events at 25 A MeV: open 
circles, full points, open squares full squares, open triangles, full triangles, open crosses refer to 
Z > 8, 18, 28, 38, 48, 58 68, respectively. 

channel is negligible. 

The close similarity between statistical models [6] and data [4,5], together with the 
isotropy of fragment emission, already suggest that these sets of data are close to a 
statistical equilibrium. 
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Figure 2. Charge distribution for central collisions (lines) compared to the QP distributions 
(symbols) in the same calorimetric energy bins. The distributions are normalized to the charge of 
the emitting systems. Left: all detected fragments, but the heaviest. Right: heaviest fragment. 
Fission has been recognized and reconstructed through a correlation technique [4, 5]. 



To progress on this point, we can also compare different data sets. Central events 
correspond to a narrow distribution of the excitation energy, while the quasi-projectile 
data supply a widely spread excitation function. Therefore it is possible to find for 
these data sets common values for the energy deposited in the source and to compare 
charge distributions and charge partitions. In the left four panels of Fig. 2 the charge 
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distribution of all the reaction products, but the largest one, masured for central Au + C 
events at 25 A MeV, Au + Cu at 25, 35 A MeV and Au + Au at 35 A MeV (top left, 
top right, bottom left, bottom right panels, respectively) are compared to the quasi- 
projectile distributions at the same excitation energy. The calorimetric values of the 
energy deposited in the systems cover a wide range. Indeed their values are 1.6, 3.1, 4.7 
and 7.4 A MeV, respectively. Moreover, the comparison among the charge distributions of 
the heaviest fragment in each event (Fig. 2) allows to check the scaling of higher moments 
like the variances, even in more detail. 

A very good scaling behaviour is also obtained by comparing the presented sets of 
data and the charge distributions measured by the FASA collaboration for the reactions 
p, a + Au at 8.1, 4 and 14.6 GeV [7] and by the Indra collaboration for Xe + Sn reaction 
at 32 A MeV incident energy [8]. 

The remarkable scaling between different sets of data means that these data samples 
can be analyzed, at least in a good first approximation, within statistical microcanonical 
methods. 

2.2. PSEUDO-CRITICAL BEHAVIOUR IN NUCLEI 

Since the early 80's, size distributions have been fitted with power laws [9], and more 
sophisticated critical analyses have been performed following theoretical concepts coming 
from percolation theory. More recently, an astonishing good scaling behaviour has been 
observed in the EOS data [10] and tentatively associated to the critical point of the 
nuclear liquid-gas phase transition expected to occur in nuclear matter in the framework 
of the Fisher droplet model [11]. The debate on the order of the transition has been 
further animated by a very recent analysis of the EOS and Isis data [12] which shows a 
high quality scaling of the fragment size distribution over a wide range of charges and 
deposited energies with an ansatz for the scaling function taken from the Fisher droplet 
model. The Fermi gas "critical" temperature (about 8 MeV) extracted in these papers 
is identified as the temperature of the thermodynamical critical point and the whole 
coexistence line of the nuclear phase diagram is reconstructed under the assumption that 
the Fisher model gives a good description of the multifragmentation phenomenon [12]. 
In this interpretation, multifragmentation would correspond to the critical point of the 
nuclear matter EoS (with a lower temperature due to finite size and Coulomb effects), i.e. 
to a second order phase transition. 

However, the experimental observation of a flattening of the caloric curve [2,13] and 
recent studies [14] in the framework of the SMM model [6] (in which the resulting critical 
temperature is about 20 MeV) point rather to a first order phase transition, and this is 
also suggested by other thermodynamical statistical multifragmentation models [15]. One 
may also wonder what physical meaning can be attributed to a Fermi gas estimator of the 
temperature at excitation energies of several A MeV. Indeed for the Au quasi-projectile 
data we are discussing, the temperature, evaluated from an isotope and a kinetic energy 
thermometer [16], strongly deviates from the Fermi gas ansatz (left part of Fig. 3). 

Surprisingly enough, a comparable quality of scaling (Fig. 3) and a consistent set of 
critical exponents can also be extracted from the Au QP data samples, for which the heat 
capacity evaluation points to a first order phase transition [4,16], and the quality of the 
scaling does not depend the ansatz adopted for the T parameter. The same ansatz for 
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Figure 3. Left side: different estimations of the QP temperature. Circles: kinetic energy 
thermometer from Eq.(3) and 3-rd degree polynomial fit (solid line). Stars: isotopic thermometer 
from Ref. [16]. Dashed line: Fermi gas ansatz y/8E*. Middle panel: QP scaled yields as a 
function of the scaled T, estimated from the solid line (squares) and from the dashed line (circles). 
Right panel: Scaled yields for central collisions. In the middle and right panels to represent the 
results on the same picture a constant horizontal shift C is applied to the distributions. 

the scaling function applied to central events [5] gives similar exponents and points to the 
same critical-like energy. 

However, critical exponents and scale invariance are compatible with many different 
physical phenomena and are not necessarily linked to a thermodynamic second order phase 
transition [17]. In particular the observed signals of critical behavior can be compatible 
with a first order phase transition, since in different statistical models size distributions 
that mimic a scale invariant behavior are observed inside the coexistence zone of small 
systems [18]. This means that scaling per se does not demonstrate the existence of a 
phase transition, and even less defines its order or allows to localize the system on the 
phase diagram. 

On the other hand, the fact that all the analyzed reactions behave as a universal 
multifragmentation process independent of the entrance channel and directly correlated 
with the available energy only, is a strong indication of a microcanonical equilibrium. 

3. MULTIFRAGMENTATION AND SIGNALS OF PHASE TRANSITION 

It has been shown [19] that for a given total energy the average partial energy stored in 
a subsystem of the microcanonical ensemble is a good thermometer while the fluctuations 
associated to the partial energy can be used to evaluate the heat capacity. An example 
of such a decomposition is given by the kinetic E k and the interaction Ej energies. In 
particular first order phase transitions are marked by singularities and negative heat 
capacities [15,19], corresponding to fluctuations anomalously larger than the canonical 
expectation. If the system is in statistical equilibrium, a measurement of anomalous 
fluctuations at a given energy is an unambiguous proof of a thermal first order phase 
transition. 

Experimentally, the total energy E* deposited in the system can be evaluated event by 
event by calorimetry. If one is also able to reasonably estimate for multifragmentation 
data the relation between the measured charge and mass of the reaction products and 
the value of the freeze-out volume, one can obtain the event-by-event Q-value and the 
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Coulomb energy, and thus the interaction energy, according to: 



M 



E I = Y,m + E F co ° ul -m Q 



1,M) 



(1) 



where M is the total multiplicity, rrii (m ) are the mass excess of the primary products 
(of the source). The interaction energy fluctuation can then be studied as a function of 
the total energy and the heat capacity can be evaluated according to: 



C 



C 2 k T 2 



C h T* 



(2) 



where a\ = erf is the fluctuation of the interaction energy Ej from Eq.(l), T is the 
temperature, and C k is the kinetic heat capacity that can be evaluated by taking the 
numerical derivative of {E k ) = E* — (Ej) with respect to T. Eq.(2) shows that a negative 
heat capacity corresponds to partial energy fluctuations in the microcanonical ensemble 
that exceed the corresponding fluctuations in the canonical ensemble (a can = CkT 2 ). 
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Figure 4. Average fragment kinetic energy for QP data. Symbols: experimental values, lines: 
many-body Coulomb trajectories for a volume 2.7 Vq (full line), 4Vb (dotted), 5Vq (dash-dotted), 
6Vb (dashed). 

The average freeze-out volume can be determined [16] from the detected fragment 
kinetic energies through a many-body Coulomb trajectory calculation (Fig. 4). For the 
sets of data here presented, a volume close to three times the normal volume reproduces 
the measured kinetic energies. Only for the central Au + Au collisions, the interplay 
between the Coulomb energy and the radial flow (about 1 A MeV) gives an uncertainty 
in the determination of the freeze-out volume. Indeed, the fragment kinetic energies would 
also be compatible with zero radial flow and an increased Coulomb repulsion from a more 
compact configuration. This ambiguity, however, does not affect our main conclusions, as 
is evident from the results shown in Fig. 5. 

An estimator of the microcanonical temperature T can be obtained by inverting the 
kinetic equation of state: 



M Ai 
( E k) = ( ~Fr\ 



)T 2 



l))T 



(3) 



The unknown parameters of Eq.(3) are the average side feeding correction A A on the 
fragment masses At, and the level density parameter a of primary fragments. Since the 
Coulomb energy is positively correlated with the charged products multiplicity, the value 



7 



obtained for the freeze-out volume in each excitation energy bin depends on the side 
feeding correction. This means that VpoiO-i an d the percentage of evaporated particles 
have to be fixed consistently with an iterative procedure. When this is done [16], we find 
that the microcanonical temperature, evaluated through Eq.(3) results in agreement with 
the measured isotope temperature (left panel of Fig. 3). 
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Figure 5. Left panel: Normalized partial energy fluctuations and kinetic heat capacity for QP 
events (grey contours) and central Au + C (black dots), Au + Cu (squares, triangles), Au + Au 
reactions before (open stars) and after subtraction of 1 A MeV radial flow (black stars). Left 
panel: Heat capacity per nucleon of the source for QP events and central reactions. 

In order to minimize spurious fluctuations due to unmeasured quantities, a constant side 
feeding correction is applied in each excitation energy bin, and no volume fluctuations 
are allowed [16]. The partial energy fluctuation overcomes the canonical expectation 
approximately at the same values of E* for peripheral and central collisions data. The 
same thing is true for central Xe + Sn collisions from 32 to 50 A MeV, measured with 
the Indra device [16]. Therefore, the different data sets are fully compatible, the results 
do not depend neither on the detector, nor on the data selection. 

It is however important to stress that the quantitative study of nuclear thermodynamics 
is still at the first stage. Only systematic studies of correlated observables, performed 
with sophisticated experimental devices, with low energy threshold, high granularity and 



S + 3S Ni (grey), + 04 Ni (hatched) 




A 

Figure 6. Mass probability distribution for a Si-Csl Garfield telescope placed at 15° in the 
reactions 32 S+ 58 Ni (grey histograms) and 32 S+ 64 Ni (hatched histograms) at 15 A MeV. 
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isotopic resolution will bring a detailed understanding of the phase transition and will 
allow to locate the position of the observed multifragmentation in the phase diagram 
of nuclear matter [1,2]. Recent investigations on the opening of the multifragmentation 
channel, on the caloric curve and on the limiting temperature for A 100 systems started 
at moderate incident energies (~ 15 A MeV) at INFN Laboratori Nazionali di Legnaro 
with the Garfield apparatus. First results [20] indicate that already at an excitation energy 
about 3 A MeV, multifragmentation events with nearly equal size fragments represent a 
non negligible fraction of the decay channels. The very good mass resolution of the 
Garfield Si-Csl telescopes (Fig.6) allows an isotope analysis. The dependence on the 
isospin of the entrance channel is apparent. Isoscaling analyses, now in progress, will 
bring information on the isospin distillation at the onset of the multifragmentation [3]. 

The authors would like to thank Garfield collaboration for kindly providing preliminary results. 
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